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I.  COORDINATE  SYSTEMS 

Fixed  at  the  center  of  the  earth  is  the  origin  of  an  XYZ  coordinate  system  (Fig.  1).  The 
Y  axis  ponses  through  the  point  at  which  the  initial  conditions  arc  established,  the  XY  plane 
contains  the  target,  and  the  Z  axis  completes  a  right-hand  coordinate  system. 


Fig.  1.  Sp^ce  Coordinate  System 


The  rocker  center  of  mass  is  locat'd  by  the  cylindrical  coordinates  r,  y,  r:  r  is  the 
distance  from  the  origin  to  the  projection  cf  the  rocket  center  of  mass  on  the  A?  plane,  y  is  the 
angle  between  Y  and  r.  and  ?  is  identical  to  Z.  There  are  two  sets  of  space  axes  located  with 
their  origins  at  the  rocket  ee:.:ct  ol  ,r.as'..  7  :»  ■_  •  '-v  -><>  oriented  that  ih:  a. 

are  parallel,  respectivelv.  to  the  XYZ  axes.  The  xyc  coordinate  systci.',  is  ss  oriented  that  <?«•' 
xy  plane  is  parallel  to  the  XY  plane,  and  y  is  parallel  to  r.  Also  located  with  its  origin  at  the 
rocket  center  of  mass  is  a  set  of  £17  £  body  axes  (big.  2).  The  roll  axis  is  cf,  the  yaw  axic.  is  r/ 
and  the  pitch  axis  is  £. 


The  orientation  of  the  body  with  respect  to  the  AgP  coordinate  system  is  described  bv 
the  Ealer  angles  4> ,  0,  if/ ,  in  the  following  manner  (Fig.  3)1.  Starting  with  C.  axes  respectively 
coincident  with  the  Agy  axes,  the  £aet  ia  rotated  counterclockwise  by  an  angle  <+>  about  the 
tfaxis.  The  second  rotation  is  counterclockwise  by  an  angle  6  abont  the  cf  axis;  whereas,  the 
third  rotation  is  counterclockwise  bv  an  angle  </>  about  the  £axis.  Also,  it  should  be  noted  that 

^Figure  4-6  from  Classical  Mechanics  by  Herbert  GoMstein,  Copyright  1950  by  Addison-Wesiev  Publishing 
Company,  Inc.,  Heading,  Mass.  t!.5. A.  tVrmisnion  to  use  this  figure  in  this  report  only  granted  by  the 
copyright  owner.  Addison-Wrsley  Publishing  Comp.ar.v,  Inc. 
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the  A/if  coordinate  system  is  obtained  by  rotating  the  xyz  coordinate  system  counterclockwise  by 
the  angle  y  about  the  2  axis. 

The  matrix  A  (y  </>,  6 ,  <fi)  is  used  to  transform  vectors  from  xy2  coordinates  to  A/if 

;  whereas,  its  inverse  A.  *  (y,  <5,  d,  tf/)  is  used  to  perform  the  inverse  transformation. 
These  matrices  are: 


A  -- 


cosV'  cos(<£  *y)  -cos#  sin(0  *y)  sin'/'  cosi/i  sint'/j  +  y)  +  cos#  cos(^  +  y)  sin/' 
-sin'/'  cos(*/>  *y)  -cos#  sin('i> +y)  cos'/'  -sin'/'  sin(d>  +-y)  +  cos#  cos(<j> +y)  Cosi/i 
sin#  sin(</>+y)  -sin#  cost'/’  +>) 


sin'/'  sin# 
0031/1  sin? 
cob# 


cos'/'  cos(</i+y)  -cos#  sin(<6+V)  sini/i 
cos'/'  sinf'/'+y)  +  ccs#  no s(0 +y)  sin'/' 
sin'/’  sin# 


-sin1/’ cos(0  4-y)  -cos#  sin(i£  +  y)  cos'/'  sin#  sin(i/>+y) 
-sinii  sin(^  ♦y)  +  cos#  cos(d>  +  y)  cosi/j  -sin#  eoaf'i +y) 
cos'/'  sin#  cos# 
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II  TRAJECTORY  MOTION 

If  motion  in  the  z  direction  iss  small,  the  Eo  '-lions  describing  the  trajectory  of  the  rocket 
center  of  mass  can  be  written  in  the  cylindrical  coordinates 


*Y 


dv%  D  vzvy 


dvy  L  vl 


~  -  So 


d  t.  S 


vx  =  component  of  velocity  vector  along  x  axis 
vy  -  component  of  velocity  vector  along  y  axis 
v .  =  component  >f  veiociiy  vector  along  z  axis 
P  component  ■>(  aerodynamic  force  vector  along  x  axis 
I*  -  component  of  aerodynamic  force  alone  v  axis 


5  •  cnrrtponer*  of  aerodynamic  force  along 
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m  =  coasting  mass  of  the  re-entry  body 
go  =  sea-level  acceleration  of  gravity 


r0  =  sea-level  radius  of  the  earth 


III.  ANGULAR  MOTION 


Euler’3  Equations  are  used  to  describe  the  angular  motion  of  the  body  about  its  center  of 


mass 


do Of 

dt 


'v - "f 

at 


Wt-V-^n 


dc,>£ 

!,  «,  </*  -  V-  *< 

dt 


>  (3) 


where 


<  j  -  principal  ruocient  of  inertia  about  the  ich  axio 
*0-  -  cyniyoBsm  of  angular  velocity  vector  a'ons-  the  ?>h  axis 


moment  along  the  tth  axis 


» <-\.o 

JWO 


t  Av 


i  *  - 


i  /  -2  a  .>-«  a 

iia  Ot  jy 
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In  order  to  obtain  the  Equations  of  motion  in  terms  of  Euler  angles.  Dew  components  I? 
tile  angular  velocity  vector  are  defined. 


In  the  body  system  of  coordinates,  cd^  has  the  components 


«  to#  siu0  sint/t  (w^,) v  -  ti ^  sin 0  cos \f>  (u^)  r  *  &>£  cond 


toy  has  t'  e  component# 


(oy)-  =  a,B  cor.>p  [t06) v  =  -  03g  sinyfr  (wy)^  =  0 


and  oj^has  the  components 


<«vV  = 0 


U’.i)  Tt  ’  0 


(coj,)  r  = 


Adding  these  componenta,  the  components  of  the  angular  velocity  vector  in  body  coordinates  are 
as  follows: 


o) ^  =  <d^  t.ia  8  sind  +  cay  cosy/ 


sin(>  cos C'  —  <dc  sin'i 


<D  •  -  CD^  COsO  -  Co^, 


After  taking  time  derivations,  it  follows  that 
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da# 

dwf 

d<»v 

sin# - 

ainy*  — — 

■-  cosip - -  a cost? 

dt 

dt 

dt 

dioe 

do)  g 

do)v 

-  = 

coslA - 

-  sinii -  —  0)^  0)^  sint? 

dt 

dt 

dt 

da ^  do  £  d 

- =  -  -  —  L<u^,  cost?] 

dt  it  dl 


Combining  Eqs.  (3)  with  Eqa.  (6)  gives 


dhi. 


aint?  — 


dt 


ntatp  cosxp 

-  [M £  +  (7^.  -/-)  CJV  <Uj]  - [Mr  +  ilr-lf)  &)r  tjf]  -  (O^ia)^  COS 9-0^) 

'i  lv 


dtiij)  cost// 

-  - [Ms  t  (7^  -/^  w,,  &.{] 

dt  I  c 


sin  ip 

- L.V^  +  (7;  - /5‘)  cj  j  ]  -  to^  a.j,  sin 0 

!v 


do ,  1  d 

-  =  —  [.V  J.  +  (/^  -  7^  )  CUc  -  -  [o^  cost?] 

dt  I,  dt 

(7) 


The  rotational  .notion  can  he  completely  described  by  Eqs.  (.4),  (5)  and  (7). 
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!Y.  AERODYNAMIC  FORCES  AND  MOMENTS 


Aerodynamic  force  and  moment  coefficients  ere  usually  available  in  body  coordinates; 
therefore,  it  is  necessary  to  find  the  components  of  the  stream  velocity  vector  in  body  coordinates 


"f 

-  (vx  -  wj 

uv 

'rr 

II 

-  vy 

' 

1 

. 

-  (l’z  -  «v 

u£  =  —  ivx  —  wx)  [cos</>  cos(<£  +  y)  —  cos0  sin(<£  +■  y)  sini^] 

-  Vy  [co80  sin(<£  +  >)  +  cos9  coe(<£  +  y)  sini/rj  -  \,v  z  -  w  z)  sini£  sin<9 

urj  *=  (nj  -  u'I)[ainy)  coa(<£  +  y)  +  cos0  sin(<5  +  y)  coin  *] 

*■  uy  [sin^  sin(<£  +■  v)  -  cos#  coa(<£  +  y)  cos  if>)  -  (i>z  -  w  t)  cost 0  sin  C 

u;  *  -  ft',  -  u>j)  sin#  sin(<£  +  y)  +  v  3in#  cos(<£  +  y)  -  (v .  -  uO  cos# 


fn  oody  coordinates,  the  aerodynamic  forces  are  then  found  from 


“(1  r 

?i  - - ?W  \U\+UV+U\\  j2  Ci('a’  5’  *'''  Re) 

u.  2  L  ■> 


»  ~  <f.  4  do) 


a.  =  angle  of  attach 
1 3  =•  roll  ancle 
3  -  fin  deflection 


Initiations  (2)  require  aerodynamic  fo’ce  components  in  space  coordinates;  therefore  F  c  .  F  and 

.  .  C  77  ’ 

£  mu«t  he  transformed  to  space  coordinates 
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D 

Ff 

L 

=  A"1 

Fv 

S 

Fl 

(11) 


Since  the  aerodynamic  coefficients  are  functions  of  angle  of  attack  and  roll  angle,  it  is  necessary 
to  determine  these  angles. 


D  =  ^  [  cos^  cosfci  +  y)  -  cos#  sin(<£  +  y)  sin^r]  -  F ^ 

[sin^  cos(<£  +  y)  +  cos#  sin(<^>  +  y)  cos 0]  +  sin#  sin(^>  +  y) 

L  ~  F^  [cos0  sin(^)  +  y)  +  cosO  cos(<£  +  y)  sin^]  —  F ^ 

[sin^  sin(^>  +  y)  -  cosfl  cos(<p  +  y)  cos0  ]  -  F  ^  sin#  cos(<£  +  y) 
S  =  F^  sin#  sin^  +  F ^  sin#  cosy)  +  F  ^  cos 0 


>  (12) 


J 


cos  a  - 


H  +  %2+  “{) 


C08  ft  = 


/..  2  ,  2  \  ^ 
V*  -7  +  u  i  j 


(13) 


The  pitching  moment  is  given  by 


"i  = 


"  %  1 


uj  2 


p(u|  +  u~  u2j  d3  CMC  (a,  ft,  8 ,  M,  /?e)  (14) 


the  yawing  moment  by 


ir  1 


W  = 

V 


/»(“|  ■*  B  +  U2AJ3  CMV  (a,  ft,  8,  k.  Re) 

O  ‘ 


“Ci  2 


(15) 
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and  the  rolling  moment  by 


(«? 


(16) 


V.  SUMMARY 

The  trajectory  Equations  (1  and  2)  simultaneously  with  the  rotation  Equations  (4,  5,  and  /) 
describe  the  dynamics  of  the  coasting  rocket.  The  aerodynamic  forces  are  given  by  Eqs.  (10)  and 
(12),  whereas,  the  aerodynamic  moments  are  obtained  from  Eqs.  (14),  (15),  and  (16).  The  forces 
and  moments  are  functions  of  angle  of  attack,  roll  angle,  and  fin  deflection.  Angle  of  attack  and 
roll  angle  (ire  computed  by  Eq.  (13);  whereas,  fin  deflection  is  a  programmed  input.  All  aerodynamic 
coefficients  must  be  supplied  as  inputs. 
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NOMENCLATURE 

transformation  matrix  from  ryz  axes  <fp£axes. 
transformation  matrix  from  £  axes  to  .ryz  axes, 
aerodynamic  coefficient  of  normal  force  in  yaw  plane, 
aerodynamic  coefficient  of  normal  force  in  pitch  plane, 
aerodynamic  coefficient  of  cord  force, 
aerodynamic  coefficient  of  pitch  about  center  of  mass, 
aerodynamic  coefficient  of  yaw  about  certer  of  mass, 
aerodynamic  coefficient  of  roll  about  center  of  mass, 
rocket  base  diameter. 

aerodynamic  force  component  along  x  axis, 
aerodynamic  normal  force  in  yaw  plane, 
aerodynamic  normal  force  in  pitch  plane, 
aerodynamic  cord  fore*', 
sea-level  acceleration  of  gravity, 
pitch  moment  of  inertia, 
yaw  moment  of  inertia, 
roll  moment  of  inertia. 

distance  between  missile  axis  and  fin  center  of  pressure, 
aerodynamic  force  component  along  y  axis. 

mass  of  coasting  rocket. 

Mach  number. 

pitch  moment,  see  Eq.  (14). 
yaw  moment,  see  Eq.  (15). 
roll  moment,  see  Eq.  (16). 
sea-level  radius  of  earth. 
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NOMENCLATURE  (Cont’d) 

Re  =  Reynolds  number. 

r,  y,  z  =  earth-fixed  secondary  coordinate  system  (Fig.  1). 

S  =  aerodynamic  force  component  along  z  axis, 
a  f ,  ,  «£  =  components  of  stream  velocity  in  body  coordinates. 

vx  =  center  of  mass  velocity  vector  component  along  x  axis. 

t>  =  center  of  mass  velocity  vector  component  along  y  axis, 

o  =  center  of  mass  velocity  vector  component  along  ;  axis. 

w  =  downrange  wind  velocity. 
w  -  crossrange  wind  velocity. 

x,  y,  z  =  local  space  coordinates  associated  with  r,  y,  z  (Fig.  1). 

X,Y  ,Z  =  earth-fixed  primary  coordinate  system  (Fig.  1). 

a  =  angle  of  attack:  angle  between  stream  velocity  vector  and  f  axis. 
ft  =  roll  angle:  angle  between  plane  of  angle  of  attack  and  rj  axis. 

5  =  fin  deflection. 

X,  11,  v  -  local  space  coordinates  associated  with  A  ,  Y  ,  Z  (Figs.  1  and  3). 
r) ,  £  -  body-fixed  coordinate  uystem  (Fig.  2). 
p(r)  =  atmospheric  density. 

<f> ,  6 ,  \f>  -  Euler  orientation  angles  (Fig.  3). 
ojj  =  pitch  rate, 
t ov  ■*  yaw  rute. 

=  roll  rate. 

"<*>  •  <*>4,  =  see  Eq.  (4). 
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